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The nature of transverse momentum spectra of the secondaries produced in various high energy nuclear and 
particle collisions!^ |[ ||, is still of much importance and interest for some various well-known reasons. 

Pions constitute the most abundant variety among the various secondaries produced in high energy collisions 
which include particle-particle, particle-nucleus and nucleus-nucleus interactions. Within the confine of the 
present work, the word ‘particle’ refers mainly to proton, because no other hadron, as either projectile or target, 
has here been reckoned with. In this work we would concentrate on analyzing the nature of transverse momentum 
dependence of the inclusive pion production cross section for some high energy nucleon- and nucleus-induced 
reactions and also on understanding the characteristics of some of the cross section-ratio versus transverse 
momentum(p 7 ’) plots. This indicates certainly no preferential treatment to pT-spectra as such; rather it is 
because of the fact that the rapidity spectra of pions have already been studied with an identical approach, 
and the encouraging results of our investigation have been communicated in a separate work|Q. The ratio- 
behaviours are significant, mainly because of the fact that they offer a dependable consistency check-up of the 
final working formula which has here been framed and forwarded mostly from the phenomenological points of 
view; and which has also finally been put into use for actual calculations of the results. Besides, comparing the 
results obtained by the model of our choice with some of the existing and the most prominent models would 
also be one of our objectives here. Still, trying to obtain first reliably good fits to the data on inclusive cross 
sections versus transverse momenta (pr) values of the pions produced in NN, NA, and AA interactions from 
a very generalized approach constitutes here our topmost concern and priority in this work. 

Regarding the very basic motivations of this work, we offer precisely the following few points: (i) The 
high energy physicists(both theorists and experimentalists) have, by and large, focused heavily on ruT-scaling 
behaviour put forward by Hagedorn for PP collision in his statistical model |8|. But for production of particles 
at relatively high transverse momenta the ruT-scaling hypothesis is seen to be experimentally inapplicable 
for which Hagedorn himself later shifted to a power law form[^ of expression which would be used in the 
present work. The universal function in the Hagedorn model with mT’-scaling is a Bose-Einstein or Fermi 
Dirac distribution and depends on a simple parameter, i.e., the slope parameter. Contrary to this picture we 
would try to attack the problem in reverse and try to determine the universal behaviour with no particular 
concern about the status of mT'-scaling by analysing intensively and extensively the experimental data on both 
the nucleon-nucleon(PP) and nucleon(nucleus)-nucleus interactions. This aspect is unique in our approach, 
(ii) Secondly, we have attempted to show somewhat clearly the degree of dependence of AA collisions on PP 
collisions and also the extent of intertwining between proton-proton and nucleus-nucleus reactions. In fact, 
the entire energy-dependence of particle production characteristics of nucleus-nucleus collisions has here been 
transmuted to the case of proton-proton reaction at the same center of mass energy. This has helped us to 
construct a pathway between NN and NA/AA interactions, (iii) Thirdly, our approach to the analysis of the 
problem is not based on any apriori assumptions of the standard and textual premises of the QGP hypotheses 
or about the constituent picture of any of the popular brand of theories. 

The organization of this paper is as follows. In the next section (Section 2) we give the outline of the 
basic outlook and the approach to be taken up for this study. In Section 3 we present the essential steps in the 
methodology of our work here. The Section 4 contains the results of model-based calculations and some broad 
discussion on the results obtained by this combinational approach. In Section 5 we have made a comparison 
of the results obtained by the model of our choice and those by some other very front-ranking ones. The last 
section is reserved for summing up the conclusions and making some points of observational interest and general 
validity. 

2. The Basic Approach : Tracing the Outline 

Following the suggestion of Faessler|l^ and the work of Peitzmann||ll| and also of Schmidt and Schukraftfl^, 
we propose here a generalized empirical relationship between the inclusive cross-section for pion production in 
nucleon(N)-nucleon(N) collision and that for nucleus(A)-nucleus(A) collision as given below: 

(AB ttX) ~ {PP ttX) , 

dp'^ dp’^ 

where px) could be expressed in the factorization form, px) = f{y) g{pt)- 


( 1 ) 
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While investigating a specific nature of dependence of the two variables {y and px) either one of these is 
assumed to remain constant. In other words, more particularly, if and when the p(-dependence is studied by 
experimental groups, the rapidity factor is treated to be constant and the vice-versa. So, the formula turns into 

^0 ^ ^ {PP - , (2) 


The main bulk of work, thus, converges to the making of an appropriate choice of form for g{pT)- And 
the necessary choices are to be made on the basis of certain premises and physical considerations which do not 
violate the canons of high energy particle interactions. 

The expression for inclusive cross-section of pions in proton-proton/antiproton scattering(see figures 1 and 
2) at high energies in Eqn.(2) could be chosen in the form suggested first by Hagedorn|P: 

^0 ^ ttX) = Cl ( 1 + , (3) 

dp'^ po 

where Ci is the normalization constant, and Po, n are interaction-dependent chosen phenomenological 
parameters for which the values are to be obtained by the method of fitting. Their -^/s-dependences are here 
proposed to be given by the following formulations: 


PoiVs) = a + 


/s In(Y^) 


( 4 ) 


where a = 1.5, b = 79.4 and 


i(a/s) = d + 


In (Ts) 


( 5 ) 


with d = 6.5, b = 127. The -^/s-dependence of po and n would be shown later diagrammatically in 
the text and their data-base would also be indicated. The nature and significance of these parameters could 
be appreciated from the work of Hagedorn|Q, and those of Bielich et al[|^ and Albrecht et al|Q. The ratio 
Po fn characterizes the intercept of the slope of the transverse momentum spectra in the limit px 0 and has 
a value near 150 MeV for PP and also for some nuclear collisions. The prescribed limit of its value is 130 — 170 
MeV which we have shown to be correct and consistent in our work. 

The final working formula for the nucleus-nucleus collisions is now being proposed here in the form given 
below: 




ttX) « C 2 {PP 

dp"^ 


ttX) 


( 6 ) 


with g{px) = (e + cipr — PPx)- This quadratic parametrization is our own suggestion and is called 
De-Bhattacharyya parametrization(DBP). In the above expression C 2 is the normalization term which has a 
dependence either on the rapidity or on the rapidity density of the pion; e, a and (3 are constants for a specific 
set of projectile and target. 


Earlier works |l4 |l5i |l6|] showed that g{px) is less than unity in the pT-domain, px < 1.5 GeV/c. Besides, 


it was also observed that the parameter e, which gives the value of g{px) at px = 0, is also less than one and this 
value differs from collision to collision. The other two parameters a and [3 essentially determine the nature of 
curvature of g{px)- However, in the present context precise determination of e is not possible for the following 
understated reasons: 

(i) To make our point let us recast the expression for (6) in the form given below: 




ttX) « C 2 {A.Bf + ^ )-" 

Po 


( 7 ) 


Quite obviously, we have adopted here the method of fitting. Now, in Eqn.(7) one finds that there are two 
constant terms C 2 and e which are neither the coefficients nor the exponent terms of any function of the 
variable, px- And as e is a constant for a specific collision at a specific energy, the product of the two terms C 2 
and {A.By appears as just a new constant. And, it is quite difficult to obtain fit-values simultaneously for two 
constants of the above types through the method of fitting. 

(ii) From Eqn.(2) the nature of g{px) can easily be determined by calculating the ratio of the logarithm 
of the ratios of nuclear-to-PP collision and the logarithm of the product AB. Thus, one can measure e from 
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the intercept of gipr) along y-axis as soon as one gets the values of for both AB collision and PP 

collision at the same c.m. energy. In the present study we have tried to consider the various collision systems 
in as many number as possible. To do so, we have to consider the data on normalized versions of for 

many collision systems for which clear if^^-data were not available to us. Furthermore, from these normalized 
versions we can/could not extract the appropriate values of E^^ as the normalization terms, total inclusive 
cross-sections((Ji„) etc., for these collision systems cannot always be readily obtained. Besides, it will also 
not be possible to get readily the data on inclusive spectra for PP collisions at all c.m.energies, like e.g., at 
-y/i = 17.8Gey(c.m. energy of Pb + Pb collision). 

In order to sidetrack these difficulties and also as an escape-route, we have concentrated almost wholly to 
the values of a and (3 for various collision systems and the effect of C 2 and e has been combined into a single 
constant term C 3 . Hence, the final expression becomes, 

E^{AB ^ ttX) « C 3 {A.B)^°‘P^-^Pt) (1 )-" ( 8 ) 

dp"^ Po 

with C 3 = C 2 {A.By. 

The exponent factor term apr — fipr obviously represents here [gipr) — e] instead of gipr) alone. Thus, 
after obtaining fit-values of a and /3, if [gipr) — e] are plotted for various collision systems, all the curves would 
originate from a single point, i.e. origin; and the systems and processes are then really comparable. In other 
words, in this convenient way we could study and check the scaling characteristics of g{pT) with respect to the 
collision systems. 

The expression( 8 ) given above is the physical embodiment of what we have termed to be the grand 
combination of models(GCM) that has been utilized here. The results of PP scattering are obtained in the 
above on the basis of eqn.(3) provided by Hagedorn’s model(HM); and the route for converting the results of 
NN to NA or AB collisions is built up by the Peitzmann’s approach(PA) represented by expression(2). The 
further input is the De-Bhattacharyya parametrization for the nature of the exponent. Thus, the GCM is the 
totality and resultant of HM, PA and the DBP, all of which are used here. 

3. The Procedure 

At the very start let us mention a few points of special concern. Though pions are produced in all three varieties 
in the physical processes, we have treated their cases here in an average and charge-independent way. Secondly, 
we treated theoretically the PP and PP reactions at accelerator or collider energies on an equivalent footing, as 
annihilation channels are found to have no or negligible contributions at high energies, > 10 GeV. Besides, 
while collecting data-sets we have, at times, assumed and used the total negative particle multiplicity to be 
physically equivalent to negative pion multiplicity for all practical purposes, as the pions are overwhelmingly 
predominant in number compared to all other particle-species. And this is done only under compulsion arising 
out of the non-availability of direct and reliable data on pion production in the relevant collision. 

The first step towards progress in the present work consists of fitting the inclusive transverse momentum 
cross sections of pions produced in both PP and PP collisions available at various high energies with the help 
of the proposed form given in Eqn.(3). The values of Pq and n at different high energies are specifically noted 
in Table 1 and Table 2 for the case of pion production alone. Next, we draw graphs to investigate closely the 
nature of -y/s-dependence of these two model-based parameters, viz, po and n. The plots of po vs. ^/s and n 
vs. i/s are presented in Fig.3. While plotting graphs we have obviously assumed that at very high energies the 
proton-proton and proton-antiproton collisions could be treated at par and with the acceptance of equivalence 
between each other. 

Hereafter, in studying the nature of pr-spectra in all nucleon-nucleus and nucleus-nucleus collisions, the 
interaction energy in all cases of nucleus-induced reactions is invariably converted first into the c.m. system 
values, that is expressed in .^snn, then the values of po and n are picked up from the already drawn graphical 
plot shown by Fig.3. While analyzing nucleus-dependence with expression( 8 ) and trying with the fit parameters 
G 3 , a and /3, we have inserted these extracted values of po and n for PP cross section term occurring in 
the expression( 8 ) given above. The values of C 3 , a, (3, y^/ndf and some other relevant information on the 
specifities for various nucleon-nucleus and nucleus-nucleus collisions are depicted in tabular form [Table-3]. The 
solid curves in all the diagrams represent the theoretical calculations. The fits, in most cases, are done within 
the pr-band 0.8 < pr < 3.0 GeV/c. And this choice is obviously attuned to the availability of the experimental 
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data in the several nuclear collisions. However, only exceptions are the cases of DD and aa interactions at 
Ehab = 480H GeV where the p-r-ranges are much wider. And this would be evident from the Fig.5. The ratio 
values are calculated simply by dividing the valid expression(8) for the corresponding collisions between any 
pair of projectile and target with some multiplicative factors shown in the several ratio-figures. No further 
complications about them are taken here into account. 


4. Model-Based Results and Discussions 


Taking the cue from the preceding section one might state that the actual start of the present work takes place 
from plotting of the po vs. ^/s graph, and of n vs. y/s graph(Fig.3), which offer us the necessary lever to apply 
the connective corridor or connector built up by Peitzmann’s approach(PA) for converting the results of nucleon- 
nucleon collision to that of nucleon-nucleus or nucleus-nucleus interaction at high energy. The total work is the 
combination of these two models and of the parametrization referred to as DBP. The nature of agreement here 
between all model-based calculations and experimental measurements is modestly fair. The reproductions of 
the data on invariant inclusive cross section vs. transverse momentum plots, right from the proton-deuteron 
collision to the heaviest lead-lead interaction, and also the data involving some Uranium-induced collisions as 
well, are presented in the diagrams from Fig.4 to Fig.9. The numerical values shown in the seventh column of 
Table 3 depict the corresponding y^/ndf values which indicate the quantitative measure of the quite satisfactory 
nature of fits in qualitative terms, save perhaps the case of sulphur-lead reaction shown in Fig.8. Let us treat 
this as a minor exception for which various reasons might be at play in the various stages of measurements. 
But the diagrams shown in Fig. 10 merit special attention. The parameters a and (3 separately demonstrate a 
steadily flat nature at values around a = 0.18 ± 0.03 and /3 = 0.035 ± 0.009 respectively with regard to the 
product of AB. This ensures a saturation for the penetrability of projectile into particle/parton-structure of the 
target nuclei. The energy(-ys)-dependence has already been inserted and there should be no double counting. 
Now follow some specific comments on RHIC(BNL)-spectra which claim, in recent times, a special status. 

The latest stir in the field of high energy particle physics is caused by the spurt of particle production 
in AuAu collisions at -xJsnn = 130 GeV in the RHIC(BNL). We have concentrated here on the inclusive 
cross section of production of neutral pions in the minimum bias events []l7| , and also for both peripheral and 
central collisions |l^. The values of a given in the fifth column of Table-3 for AuAu collisions do not show 
any appreciable change in magnitude for the three separate cases of the event-sets, while those for (3 remain 
same within errors. So, within the limits of 0.8 < pt ^ 3.0 GeV/c, i.e. within the clear domain of ‘soft’ 
particle production the results on AuAu collisions at RHIC experiments modestly corroborate our claim that 
the parameters, a and (3 are independent of the collision centrality. However, with growing values of the 
transverse momentum(p 7 ’ > 3 GeV/c) i.e., for hard collision, the behaviour, in all probability, is likely to turn 
more complex; and initial indications to such complications are somewhat manifest by the slight but detectable 
deviations of the fits provided by us here for values around pt > 3.5 GeV/c. So, the present study does neither 
essentially contradict the content of Adcox et al]^, nor could essentially support the signal message received 
from them. Indeed, we have left aside here the measurement on average productions of charged hadron which 
obviously include other varieties of charged mesons and baryon-antibaryons as well. As their incorporations 
might spoil the exclusivity of pion production case, we have purposefully sidetracked those data in the present 
work; and we have converged more particularly to the transverse momentum spectra of only the neutral pions. 

The features of the cross section ratios(CSR) have been displayed in the diagrams from Fig.11 to Fig.16. 
The plots presented in Fig.16 on prediction alone with no data are the exceptional ones. The patterns of the 
plots of the specific ratios versus the transverse momentum ranges are, in most cases, modestly similar with 
exceptions for the cases of -p^(Fig.l2),-g^(Fig.l2) and ^^(Fig.l4) wherein the curvatures of both the data 
plots and the theoretical curves are just the opposite. But, in so far as the question of agreement between 
the theoretical plot versus data is concerned, there is no problem. The GSR behaviours in the two particular 
cases of (Fig.12), (Fig.l6a) and (Fig. 16b) are also of obverse type, though for the former the 

nature of fit between theory and experiment is quite satisfactory; and for the latter two this point does not 
just arise as they are predictive theoretical plots. But, the question of similarity or dissimilarity between any 
two or among the graphs cannot be delved into in any further detail, as the data-points shown in the graphs 
indicate the measurements at very different energies. The actual data on predicted ratio behaviour in ^pjjjp^ 
could be obtained separately from the RHIC in future. The nature of the ratio of could be of interest 

in the sense that the deuterons are the lightest known nuclei used in colliding systems and Pb is the available 
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heaviest nuclear species which is nowadays being frequently used as both projectile and target in the CERN 
SPS and will be used in future in the LHC. The dotted curves in both Fig.16(a) and Fig.16(b) indicate the 
range of theoretical uncertainties due to errors of averaged a and (3. 

The fair agreement between the data and calculations for the ratios offers us a modestly reliable and 
somewhat convincing cross-check of the model for pion production in nuclear collisions that has here been 
proposed and developed to a certain extent. 

Now, we hold the view that our analysis and results could be made far more rigorous, should the data 
on pion production in both PP and N A/A A interactions for different pT in exactly the same rapidity region 
would be simultaneously available. But by all indications, such expectations at this stage are too unrealistic 
for which we have had to remain satisfied with whatever is at hand now. It is very encouraging to note that 
the fits reproduce data on pion production even in the not-very compatible region of rapidity in PP and AA 
collisions. So, it is only natural and reasonable to summarize that the same or compatible region of rapidity for 
both types of collisions could offer only better and much-improved fits. 


5. The pT-Spectra in Nuclear Collision: A Brief Comparison of Some Models 


Based on the analysis of the experimental data, let us mention the following few points as the summary of 
the main and basic features of the pr-spectra, especially of large-pT variety, of the secondaries produced in 
nuclear collisions: (i) strong dependence on both target and projectile mass or mass number, (ii) very weak or 
no dependence on the rapidity, (iii) very weak or no dependence on the c.m.energy, (iv) dependence on impact 
parameter for those restricted collision wherein the full overlap of target and projectile does not occur, (v) no 
further Et dependence beyond the above point. With considerations on some of these characteristics, the NA34 
collaboration|l^ demonstrated that the target and projectile dependence could be treated in a similar manner 
when parametrized by a power law of the form: 


^(P + A) = ^(P + P) 

dpT dpT 

^{B + A) = ^(P + A) 

dpT dpT 


(9) 

( 10 ) 


Schmidt and Schukraft[p^ pointed out how a fit of ^{pr) to the original Fermilab P -|- A data could lead to 
the reasonably nice description of CIpt) in a two-step manner. The similarity between ^{pt) and CIpt) was 
emphasized and was interpreted as dependence of the Cronin effect on the projectile mass number analogous 
to the one on the target. 

Next, the WA80 Collaborationmade attempt to unify functionally this two-stage transitions into a 
single-step affair in the form of the eqn.(2) dealt with previously in Section 2. Therein the expression for 
9 {pt) will have to be chosen phenomenologically in such an appropriate manner as to fit the data as nicely as 
possible. The WA80 team proposed a simple linear form of fit for g^pr)- in the present work we have taken a 
polynomial(quadratic) form instead. There is one more major difference between these two approaches. The 
WA80 team made choices of po and n values for PP and AB reactions completely separately. On the contrary 
we have introduced the results of po and n directly from our studies PP/PP reactions at various energies. In 
the approach presented and pursued here the point is to find out the energy of interactions between the nuclei; 
and once it is known, the po and n values of the PP reactions at this specific energy are to be ascertained from 
Fig.3 and inducted in the calculations. 

The performance by WA80 Collaboration|l^ is essentially grounded in a sort of the physics of 
purturbative quantum chromodynamics(PQCD) approach in Hagedorn’s form, and is also embedded in string 
and hydrodynamic considerations. This combination is called here the Mixed Model(MM). So, we represent 
the calculated results of WA80 collaboration as the product of a Mixed Model(MM). The comparison of GCM 
and MM with experimental data on SS and SAu reactions is given in Fig.17a. The solid lines in all the figures 
describe results obtained by GCM. Over the entire pT-range the agreement with this MM is quite satisfactory; 
and at very high p^-values the results based on MM have shown much better agreement than those obtained with 
the GCM. But there is a change in the scenario in cases of data interpretation on PbPb and PbNb cases(Fig.l7a). 
The comparison, in this two cases, is with calculations made by Thermal Model|^ which does not show quite 
fair agreement with data for the region beyond px ~ 2.5 GeV/c, whereas the GCM addresses competently all 
the measured data-points. In Fig.17(b) a comparison has been made between two sets of models on data from 
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AuAu collisions with two different detectors, for PbPb interaction and also for OAu reaction. And the sets of 
the models are^ Modified Hydrodynamic model of Peitzmann^^ and the GCM and (ii) the PQCD-based 
HIJING model|^ vs. the GGM. The agreement with both data and the most popular versions of the two 
models could, on the whole, be rated to be modestly fair. The HIJING model is old and well-known among the 
Monte-Garlo simulation-based models. And, unlike the case of hydrodynamic model, there has not been any 
new updating of it in the recent past. 

Gompared to the GGM, the version of hydrodynamic model offered very recently by Peitzmann[^ provides 
certainly a little better agreement, especially at large values of pr- This is observed in Fig.17b. The data 
measured from AuAu reactions with two different detectors and other accessories have been shown separately 
and the other set of data is from PbPb collision. But the hydrodynamic model of Peitzmann is based on a version 
of this model of Weidemann and Heinz who includes transverse flow and resonance decays. Besides, Peitzmann 
also introduced the spatial distribution of Woods-Saxon type instead of simple Gaussian; reckoned with the 
difference between chemical freeze-out and kinetic freeze-out temperatures instead of one universal freeze-out 
temperature, and also the role of baryonic chemical potential. The certain degree of rigour in calculations is 
reflected in fair achievements in the collisions of the two heavies. But even with such physics considerations 
for actual calculations, the number of hand-inserted parameters are too large in the Peitzmann’s latest model 
compared to just three in total in our approach. In our case the entire collective effects arising out of multiple 
nucleon-nucleon collisions and the effects of nuclear geometry are absorbed by the very simple term, {AB)^^P'^\ 
And the rest depends on just the results of PP collisions. So, our approach is much straightforward compared 
to the existing other numerous complicated approaches of the models. 

6. Summary and Outlook 

What we have done here seems, by now, to be pretty clear. We have precisely introduced a simple 
parametrization of pion transverse momentum spectra for nucleus-nucleus collisions at high energies in terms of 
that for proton-proton collision at the same energy. Analyzing a large variety of colliding systems at different 
energies we have shown that the ratio of the two transverse momentum spectra approximately proportional to 
the product of the mass numbers of the colliding systems raised to a power that is a quadratic function of the 
transverse momentum of the particle under consideration. Furthermore, this quadratic function is essentially 
universal, i.e. its coefficients have only very weak dependence on the colliding systems and energies. 

Through the massive demonstration of the fair agreements between measured data and our calculations the 
merit of the totality of the present approach is modestly well-made from the functional point of view. Despite 
this, let us sum up first here the important finer physical points which are the distinctive traits of the GCM 
alone, (i) Firstly, the GCM provides the most economical description of the vast body of experimental data on 
pion production in diverse set of reactions with only three arbitrary parameters, (ii) Secondly, the choice of the 
quadratic parametrizations for pr-dependence of nuclear cross-sections is one of our important contributions, 
(iii) Thirdly, we have suggested independently an approach for computing the values of inclusive cross sections 
for production of pions in PP collision with the help of the plots on pq vs. y/s and n vs. -^/s at any high energy 
point of the c.m.energy values, (iv) Once the energy-dependence in PP collision is taken care of by this method, 
the inclusive cross sections for pion production in NA/AB interactions becomes essentially energy-independent. 
In other words, energy-dependence in nuclear collisions at high energy is manifested only through the very basic 
interacting sets of proton-proton collisions. No reflection of energy-dependence is detected in this mechanism in 
the nuclear geometry of collisions. This is in striking contrast with the most of the existing models, (v) We have 
also arrived at a sort of mass-number scaling of the parameters in the exponent of the nuclear dependent terms. 
This is also in sharp variance with the standard frameworks, though admittedly such a scenario is possible 
with only suitable adjustment of the normalization terms in the inclusive cross sections, (vi) Besides, we have 
tried to provide substantial support to the proposed connector here between NN and NA/AB reactions by 
making our studies as extensive as possible in order to impart a good degree of reliability to this combinational 
approach, (vii) Finally, this work brings out the features of ‘universality’ of all high energy particle and nuclear 
interactions in a satisfactory manner. 

But the approach, despite its successes, suffers from the following gross limitations: (i) The master formula 
that is at the root of the success is not and cannot be derived right now from any of the first principles of basic 
physics, (ii) Data on production of a specific variety of particles in PP collisions at a few high energies are an 
indispensable necessity for analyzing and explaining the case of production of this particular particle species in 
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any nuclear collision. This PP-dependency is, by all indications, an inseparable part of this grand combination 
of models (GCM), though this might appear, at times, to be an intractable problem arising out of the lack 
of measured data on PP reaction at any particular energy, (iii) Thirdly and finally, the approach is totally 
insensitive to the physics of space-time-evolution of the collisions. But, seen from the angle of calculations on 
production of the final particles, this might be considered an advantage as well. 
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i) PP Collision: 


Table 1: Fit Values of po and n for PP Collisions at different energies. 



Relevant collision-specifics 

Cl 

PO 

n 

ndf 

poln 

23 


192 ± 6 

2.48 ±0.12 

18.88 ±0.76 

1.539 

0.131 

31 

7T~ , 

219 ± 5 

2.229 ± 0.004 

17.08 ±0.04 

1.157 

0.130 

45 

ycm — 0, 

241 ± 5 

1.89 ±0.01 

14.97 ± 0.06 

1.352 

0.126 

53 

min. bias 

240 ±7 

1.83 ±0.01 

14.25 ±0.03 

2.72 

0.128 

63 


248 ±6 

1.78 ± .01 

14.07 ± 0.02 

1.621 

0.127 


ii) PP Collision: 


Table 2: Fit Values of po and n for PP Collisions at different energies. 



Relevant collision-specifics 

Cl 

PO 

n 

X^ 

ndf 

poln 

200 

charged, r}cm = 0, min. bias 

414 ±7 

1.52 ±0.01 

11.18 ±0.02 

1.614 

0.136 

500 

charged, r)cm = 0, min. bias 

327 ± 7 

1.51 ± 0.05 

9.84 ± 0.03 

1.787 

0.153 

630 

{h^ + /i“), —3.0 < rjcm < 3.0, min. bias 

302 ± 5 

1.49 ±0.01 

9.41 ± 0.01 

1.694 

0.158 

900 

charged, rjcm = 0, min. bias 

412 ± 9 

1.50 ±0.05 

9.43 ±0.02 

1.61 

0.159 

1800 

charged, —1.0 < ycm < 1-0, min. bias 

370 ± 7 

1.48 ±0.01 

8.79 ± 0.03 

1.694 

0.168 
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Table 3: Numerical Values of the parameters: a and l3 for different high energy collisions. 


Collision 

E(GeV) 

Relevant 

C3 

a 

13 

X^/ndf 

Fig. 



collision-specifics 


(GeV/c)-i 

(GeV/c)-i 


No. 

P + D 

400 

7r~, min. bias 

400 ± 20 

0.21 ± 0.06 

0.04 ±0.01 

0.005 

4 

P + Be 

400 

7r~, min. bias 

1700 ± 80 

0.20 ± 0.05 

0.040 ±0.015 

0.008 

4 

P + Ti 

400 

7r~, min. bias 

(4.5 ±0.5) X 103 

0.21 ± 0.05 

0.04 ±0.015 

0.002 

4 

P + W 

200 

7r~, min. bias 

(3.8 ±0.5) X 10^ 

0.20 ± 0.04 

0.04 ± 0.02 

1.624 

4 

P + W 

300 

7r“, min. bias 

(4.7 ± 0.3) X 10^ 

0.17 ± 0.02 

0.031 ± 0.002 

1.15 

4 

P + W 

400 

7r“, min. bias 

(9.4 ±0.6) X 10^ 

0.14 ± 0.03 

0.02 ±0.01 

0.013 

4 

P + Au 

200 

TT^, 1.5 < ^ < 2.1, min. bias 

(1.7 ± 0.2) X 10^ 

0.16 ± 0.01 

0.04 ±0.01 

0.570 

4 

P + U 

200 

(7r+ + K+), 3.0 < y < 3.9 

(3.8 ±0.1) X 10^ 

0.20 ± 0.02 

0.032 ± 0.004 

0.903 

4 

D + D 

480A 

\ycm\ < 0.4 

360 ± 16 

0.25 ± 0.05 

0.026 ± 0.008 

1.438 

5 

D + D 

1400A 

TT^ 

(24 ± 4) X 10"“ 

0.18 ± 0.03 

0.04 ±0.01 

1.627 

5 

D + D 

1900A 

TT^ 

(28 ± 4^ 10"* 

0.15 ± 0.02 

0.04 ±0.01 

1.627 

5 

D + Au 

200A 

h~, 2.0 <y < 3.0, central 

138 ± 20 

0.23 ± 0.02 

0.04 ± 0.02 

1.720 

8 

a + a 

480A 

\ycm\ < 0.4 

624 ± 50 

0.27 ± 0.02 

0.026 ± 0.002 

0.015 

5 

O + C 

200A 

TT^, 1.5 < 77 < 2.1, min. bias 

(7± 1) X 103 

0.13 ± 0.02 

0.07 ± 0.02 

1.475 

6 

o + w 

200A 

negative particles, 

1.0 < y < 1.9, min. bias 

(1.3 ±0.4) X 10^ 

0.18 ± 0.02 

0.031 ± 0.005 

1.220 

6 

O + Au 

200A 

TT^, 1.5 < 77 < 2.1, min. bias 

(9± 1) X 10"^ 

0.161 ± 0.004 

0.038 ± 0.003 

0.030 

6 

O + Au 

200A 

h~, 2.0 <y < 3.0, central 

(1.2 ±0.2) X 10^ 

0.15 ± 0.02 

0.040 ± 0.005 

1.562 

8 

O + U 

200A 

(7r+ + K+), 3.0 <y < 3.9 

(5.3 ±0.2) X 10® 

0.180 ±0.002 

0.036 ± 0.002 

1.359 

6 

s + s 

200A 

tt'^, 2.1 < y < 2.9, min. bias 

(6.1 ± 0.4) X 10^ 

0.18 ± 0.01 

0.046 ± 0.002 

1.603 

7 

S + Ag 

200A 

h~, 2.0 <y < 3.0, central 

(1.4 ±0.1) X 10® 

0.150 ±0.004 

0.04 ±0.01 

1.790 

8 

S + W 

200A 

negative particles, 

1.0 < y < 1.9, min. bias 

(1.13 ± 0.08) X 10® 

0.17 ± 0.03 

0.031 ± 0.005 

1.566 

7 

S + Au 

200A 

tt'^, 2.1 < y < 2.9, min. bias 

(2.0 ±0.2) X 10® 

0.19 ± 0.02 

0.039 ± 0.005 

1.080 

7 

S + Pb 

200A 

negative particles, 

2 < y < 4, central 

230 ± 40 

0.15 ± 0.05 

0.05 ± 0.02 

2.646 

8 

s + u 

200A 

(7r+ + _fC+), 3.0 <y < 3.9 

(9.0 ±0.5) X 10® 

0.17 ± 0.02 

0.03 ±0.01 

0.433 

7 

Au + Au 

8450A 

TT^, \ri\ < 0.35, min. bias 

120 ± 40 

0.15 ± 0.02 

0.036 ± 0.003 

0.449 

9 

Au + Au 

8450A 

tt'^, \ri\ < 0.35, central 

500 ± 30 

0.145 ±0.008 

0.036 ± 0.002 

1.524 

9 

Au + Au 

8450A 

TT^, \ri\ < 0.35, peripheral 

27 ± 7 

0.14 ± 0.02 

0.036 ± 0.002 

1.836 

9 

Pb + Nb 

160A 

TT^, 2.3 < y < 3.0, min. bias 

(2.7 ± 0.4) X 10® 

0.17 ± 0.02 

0.036 ± 0.004 

1.614 

7 

Pb +Au 

160A 

charged pions, 

2.1 < 7] < 2.6, central 

1230 ± 5 

0.18 ± 0.03 

0.035 ± 0.002 

1.822 

8 

Pb + Pb 

160A 

tt'^, 2.3 < y < 3.0, min. bias 

(9± 1) X 10® 

0.16 ± 0.03 

0.03 ±0.01 

1.070 

7 
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,3 

Figure 1: Plot of vs. pr for secondary pions produced in P + P collisions at different c.m. energies. The 
various experimental points are taken from Ref.|^^ and Ref. |2^. The solid curves give the theoretical fits on 
the basis of Hagedorn’s model(Eqn.3). 



Figure 2: The inclusive spectra for secondary pions produced in P + P collisions at y/s = 200, 500, 630,900 and 
1800 GeV. The various experimental points are from Ref.[^, Ref.[^ and Ref, . The solid curvilinear lines 
are drawn on the basis of Eqn.3. 
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Figure 3: Values of po and n as a function of c.m. energy y/s. Various data points are taken from Table-1 and 
Table-2. The solid curves are drawn on the basis of Eqn.4 and Eqn.5. 

.25 
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Figure 4: Plot of for production of secondary pions in some Proton induced collisions at E^ab = 200, 300 
and 400 GeV as a function of transverse momentum pT- The various experimental data-points are taken from 
Ref. Ip^ , Ref.p^ and Ref.[^. The solid curves represent the fits on the basis of Eqn.8. 
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,3 

Figure 5: Plot of vs. pt for secondary pion produced in D D and a a collisions at different energies. 
The various experimental points are taken from Ref.j^] and Ref. . The solid curves provide the fits on the 
basis of Eqn.8. 
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Figure 6: Nature of inclusive spectra for production of secondary pions in different Oxygen induced collisions 
at pLab = 200A GeV. The various experimental points are taken from Ref. , Ref.|^ and Ref. H ■ The solid 
curvilinear lines represent the fits on the basis of Eqn.8. 
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Figure 7: Plot of inclusive spectra vs. pt for production of secondary pions in some sulphur-induced and lead- 
induced reactions at E^ab = 200A GeV and at E^ab = 160A. The various experimental data-points are from 
Ref.[0, Ref., Ref.||^ and Ref.|^. The solid curves provide the fits on the basis of Eqn.8. 
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Figure 8: The nature of for production of secondary pions in D-\- Au, O -h Au, S-\-Ag and SPb collisions 

CiPrp 

at E^ab = 200A GeV and in Pb + Au collision at E^ab = 160A GeV with change of transverse momentum. The 
various experimental points are taken from Ref. , Ref. [Q and |Q . The solid curves represent the fits on the 
basis of Eqn.8. 
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Figure 9: Invariant spectra for secondary pions produced in Au + Au reaction at ^snn = 130 GeV(RHIC). The 
various experimental points are from Ref.and Ref. |^. Against the general background of pion production 
in Au + Au reaction, the symbols PbSc and PbGl indicate simply the two separate calorimeters; the former 
is the presentations of data-set obtained by lead-scintillator sampling calorimeter and the latter are those by 
lead-glass Cerenkov calorimeter. The solid curves give the fits on the basis of Eqn.8. 



A.B 

Figure 10: Values of a and (3 for different collisions as functions of the product of mass numbers(Ai?) of the 
interacting nuclei. The fitted values of a and /3, enlisted in TableS, are taken as the data points; and are 
denoted by empty and filled squares respectively. The dashed line gives the average values 0.18 ±0.03 for a and 
0.035 ± 0.009 for f3. 
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Figure 11: Plot of ratios of cross-sections for different Proton-induced collisions at different energies to that for 
Proton-Proton collisions at the corresponding energies. The data-points are taken from Ref. |l^. The present 
model-based results are shown by the solid curves. 



Figure 12: Plot of some cross section ratios as a function of pt- The data points are taken from Ref.[^. The 
solid curves depict the present model-based results. 
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Figure 13: The nature of cross section ratios for and p^^ with respect to transverse momentum pT- The 
data points are taken from Ref. [^, The present model-based results are shown by the solid curvilinear 
lines. 
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Figure 14: Plot of p^p , and cross section ratios with respect to transverse momentum pr- The 
data-points are taken from Ref.|^. The solid curves provide the present model-based results. 
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Figure 15: Plot of ratios of cross-sections for different Nucleus-Nucleus collisions at different energies to that for 
Proton-Proton collisions at the corresponding energies. The data points are taken from Ref.[Q ^ |^, |3^. To 
extract the ratio-data, we have divided the data for Pb + Pb collision at 160^ GeV by those for P -|- P 

collision at 200A GeV. And this is an approximation. The model-based results of the present work are shown 
by the solid curves. 


(S+S)/(P+P) ratio[ELai,=200A GeV]x1 O'® h-BH j 
- (0+Au)/(P+P) ratio[ELai,=200AGeV]x10'® ^ - 
; (S+Au)/(P+P) ratio[ELai,=200AGeV]x10''' ^ ; 
r (Pb+Pb)/(P+P) ratio[ELat,=160AGeV]x10° ^ -j 

GCM — : 


f 





Figure 16: Predictive nature of and cross section ratio with respect to the transverse 

momentum(p 7 ’) on the basis of present model. The values of the ratio are certainly not exact, unless the 
e-value for pion production in corresponding reactions at any definite energy could be ascertained. The dotted 
curves depict the uncertainties arising out of the errors in a and [3. 
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Figure 17: Comparison of results obtained on the basis of various theoretical models with those attained by 
GCM against measured data on pion production in some selected nucleus-nucleus collisions. For Fig.(17a) the 
data sources for SS and SAu, for PbPh and PbNb a re ]l^ an d jl^ respectively. The theoretical plots based 
on Mixed Model and Thermal model are found from ]14| and |16|. In Fig. 17(b) the data on AuAu and PbPb 
are obtained from [|^ and . The source of the reports on measurements for OAu reaction is Ref. |Q. The 
theoretical plots based on MHM and on PQCD-HIJING model are assorted form and [pTf. 








